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We report on the lattice quark propagator using standard and improved Staggered quark actions, 
with the standard, Wilson gauge action. The standard Kogut-Susskind action has errors of 0(a 2 ) 
while the "Asqtad" action has 0(a 4 ), 0(a 2 g 2 ) errors. The quark propagator is interesting for 
studying the phenomenon of dynamical chiral symmetry breaking and as a test-bed for improvement. 
Gauge dependent quantities from lattice simulations may be affected by Gribov copies. We explore 
this by studying the quark propagator in both Landau and Laplacian gauges. Landau and Laplacian 
gauges are found to produce very similar results for the quark propagator. 



I. INTRODUCTION 



II. LATTICE QUARK PROPAGATOR 
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The quark propagator lies at the heart of most QCD 
physics. In the low momentum region it exhibits dy- 
namical chiral symmetry breaking (which cannot be seen 
from perturbation theory) and at high momentum can 
be used to extract the running quark mass jj], || (which 
cannot be extracted directly from experiment). In lat- 
tice QCD, quark propagators are tied together to extract 
hadron masses. Lattice gauge theory provides a way to 
study the quark propagator nonperturbatively, possibly 
as a way of calculating the chiral condensate and Aqcd , 
and in turn, such a study can provide technical insight 
into lattice gauge theory. 

We study the quark propagator using the Kogut- 
Susskind (KS) fermion action, which has 0(a 2 ) errors, 
and an improved staggered action, Asqtad [0, which 
has errors of C(a 4 ), 0(a 2 g 2 ). These choices comple- 
ment other studies using Clover jl| || and Overlap || 
quarks. We are required to gauge fix and we choose the 
ever popular Landau gauge and the interesting Laplacian 
gauge [M Laplacian gauge fixing is an unambiguous 
gauge fixing and, although it is difficult to understand 
perturbatively, it is equivalent to Landau gauge in the 
asymptotic region. It has been used to study the gluon 
propagator §00. 

In SU(N) there are various ways to implement a Lapla- 
cian gauge fixing. Three varieties of Laplacian gauge fix- 
ing are used, and these form three different, but related 
gauges. This is briefly discussed in section HI. For a 
more detailed discussion, see Rcf. [O. 



In the continuum, Lorentz invariance allows us to de- 
compose the full propagator into Dirac vector and scalar 
pieces 

S- 1 (p 2 )=iA(p 2 ) 1 .p + B(p*) (I) 

or, alternatively, 

S- 1 (p 2 ) = Z- 1 (p 2 )[ ll -p + M(p 2 )}. (2) 

This is the bare propagator which, once regularized, is re- 
lated to the renormalized propagator through the renor- 
malization constant 



S(a;p 2 ) = Z 2 (a;n)S ien (n;p 2 ), 



(3) 



where a is some regularization parameter, e.g. lattice 
spacing. Asymptotic freedom implies that, as p 2 — > oo, 
S(p 2 ) reduces to the free propagator 

S^ip 2 ) ->z 7 -p + m , (4) 

where too is the bare quark mass. 

The tadpole improved, tree-level form of the KS quark 
propagator is 

s al Cp; m ) = u oi (%) Q/ 3 sin W + ™W ( 5 ) 



where p M is the discrete lattice momentum given by 
27m,, 



aL,, 



(6) 



For the tadpole factor, we employ the plaquette measure, 



The quark propagator was calculated on 80, I6 3 x 32 
configurations generated with the standard Wilson gluon 
action at (3 = 5.85 (a = 0.130 fm)|]T^|. We have used six 
quark masses: am = 0.075, 0.0625, 0.05, 0.0375, 0.025 
and 0.0125 (114 to 19 MeV). 



u = (iReTr(P^) 



(7) 



We give a detailed discussion of the notation used for the 
staggered quark actions in an appendix. As a convienient 
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short-hand we define a new momentum variable for the 
KS quark propagator, 

<7 M = sin(p M ). (8) 

We can then decompose the inverse propagator |1S| 



M(q) 



1 



16N c iq 2 
1 



— Trtf -qS- 1 } 



16N, 



Tr^ 1 }, 



(9) 
(10) 



where the factor of 16 comes from the trace over the 
spin-flavor indices of the staggered quarks and N c from 
the trace over color. Comparing Eqs. (||) and (||) we see 
that dividing out q 2 in Eq. (^) is analagous to dividing 
out p 2 in the continuum and ensures that that Z has 
the correct asymptotic behavior. So by considering the 
propagator as a function of q^, we ensure that the lattice 
quark propagator has the correct tree- level form, i.e., 



1 



i^f ■ q + m 



(11) 



and hopefully better approximates its continuum behav- 
ior. This is the same philosophy that has been used in 
studies of the gluon propagator JlJJ and Eq. (Js[) was used 
to define the momentum in Ref. [E|. 

The Asqtad quark action Q is a fat-link Staggered 
action using three-link, five-link and seven-link staples 
to minimize flavor changing interactions along with the 
three-link Naik term [jl3| (to correct the dispersion rela- 
tion) and planar five- link Lepage term |l4|] (to correct the 
IR). The coefficients are tadpole improved and tuned to 
remove all tree-level 0(a 2 ) errors. This action was moti- 
vated by the desire to improve flavor symmetry, but has 
also been reported to have good rotational properties. 

The quark propagator with this action has the tree- 
level form 

^(P! TO ) = u o« (%) a p sin (^) [l+g sin2 (^)] +m~6 a p, 
n 

(12) 

so we repeat the above analysis, this time defining 



q^ = sin(p M ) [l + g sin 2 (p M )] . 



(13) 



Finally, it should be noted that both actions get contri- 
butions from tadpoles, which can be seen in the tree-level 
behaviors of the two invariants 
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(14) 
(15) 



so inserting the tadpole factors provides the correct nor- 
malization. 



III. GAUGE FIXING 

We consider the quark propagator in Landau and 
Laplacian gauges. Landau gauge fixing is performed by 
enforcing the Lorentz gauge condition, ^^d^A^x) = 
on a configuration by configuration basis. This is 
achieved by maximizing the functional, 



F=^Tr{U^x) + Ul(x)}, 



(16) 



X . fx 



by, in this case, a Fourier accelerated, steepest-descents 
algorithm [ fl5| . There are, in general, many such maxima 
and these are called lattice Gribov copies. While this 
ambiguity has produced no identified artefacts in QCD, 
in principle it remains a source of uncontrolled systematic 
error. 

Laplacian gauge is a nonlinear gauge fixing that re- 
spects rotational invariance, has been seen to be smooth, 
yet is free of Gribov ambiguity. It is also computation- 
ally cheaper then Landau gauge. There is, however, more 
than one way of obtaining such a gauge fixing in SU (N) . 
The three implementations of Laplacian gauge fixing em- 
ployed here are (in our notation): 

1. d 2 (l) gauge (QR decomposition), used by Alexan- 
drou et al. KM. 

2. c? 2 (II) gauge, where the Laplacian gauge transfor- 
mation is projected onto SU (3) by maximising its 
trace @. 

3. c? 2 (Ill) gauge (Polar decomposition), the original 
prescription described in Ref. Q] and tested in 
Ref. §. 

All three versions reduce to the same gauge in SU(2). For 
a more detailed discussion, see Ref. 101. 



IV. ANALYSIS OF LATTICE ARTEFACTS 
A. Tree-level Correction 

As mentioned above, the idea of "kinematic" or "tree- 
level" correction has been used widely in studies of the 
gluon propagator Q and the quark propagator [Q, [|, ||] 
and we investigate its application to our quark propa- 
gators. For the moment we shall restrict ourselves to 
Landau gauge. To help us understand the lattice arte- 
facts, we separate the data into momenta lying entirely 
on a spatial cartesian direction (squares), along the tem- 
poral direction (triangles), the four-diagonal (diamonds) 
or some other combination of directions (circles). 

The Z function is plotted for the KS action in Fig. [IL 
comparing the results using p and q. In the top of Fig . yj 
we see substantial hypercubic artefacts (in particular look 
at the difference between the diamond and the triangle at 
around 2.5 GeV). We can suggest that this is caused by 
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FIG. 1: Z function for quark mass ma = 0.05 (m ~ 76 MeV) 
for the KS action in Landau gauge. Top figure is plotted using 
the standard lattice momentum, p and the bottom uses the 
"action" momentum, q. Note that this choice affects only the 
horizontal scale. 



FIG. 2: Z function for quark mass ma = 0.05 (m ~ 76 MeV) 
for the Asqtad action in Landau gauge. Top figure is plotted 
using the standard lattice momentum, p and the bottom uses 
the "action" momentum, q. 



violation of rotational symmetry because the agreement 
between triangles and squares suggests that finite volume 
effects are small in the region of interest. In the plot 
below, where q has been used, we see some restoration of 
rotational symmetry. 

The same study is made for the Asqtad action in Fig. |^. 
In both cases this action shows a substantial improve- 
ment over the KS action, and when we plot using the 
momentum defined by the action, rotational asymmetry 
is reduced to the level of the statistical errors. 

It is less clear which momentum variable should be 
used for the mass function, so for consistency we use q, 
as for the Z function. The effect of this is shown in Fig. [|. 
For ease of comparison, both sets of data have been cylin- 
der cut jl2|]. In the case of the mass function, the choice 
of momentum will actually make little difference to our 
results. 



B. Comparison of the actions 

In Fig. |]the mass function is plotted, in Landau gauge, 
for both actions with quark mass ma = 0.05. This time 
there have been no data cuts. We see that the KS action 
gives a much larger value for M(0) than the Asqtad action 
and is slower to approach asymptotic behavior. Asqtad 
also shows slightly better rotational symmetry. 

Looking back at Figs. [I] and ^ we see that the Asq- 
tad action displays clearly better rotational symmetry in 
the quark Z function and, curiously, improved infrared 
behavior as well. The Asqtad action also displays a bet- 
ter approach to asymptotic behavior, approaching one 
in the ultraviolet. The relative improvement increases 
as the quark mass decreases. In Fig. || we compare the 
mass function for the two actions at ma — 0.0125, the 
lowest mass studied here. The low quark mass has in- 
troduced less noise into the propagator with the Asqtad 
action than with the KS action. 
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FIG. 3: The quark mass function for quark mass ma = 0.05 
(m ~ 76 MeV) for the KS action (open circles) and the Asqtad 
action (solid triangles) in Landau gauge. Top figure is plotted 
using the standard lattice momentum, and the bottom using 
the "action" momentum. 



C. Comparitive performance of Landau and 
Laplacian gauges 

Fig. ^| shows the mass function for the Asqtad action 
in 9 2 (I) and <9 2 (II) gauges and it should be compared 
with the equivalent Landau gauge result in Fig. |. We 
see firstly that these three gauges give very similar re- 
sults (we shall investigate this in more detail later) and 
secondly that they give similar performance in terms of 
rotational symmetry and statistical noise. Looking more 
closely, we can see that the Landau gauge gives a slightly 
cleaner signal at this lattice spacing. 

Landau gauge seems to respond somewhat better than 
c? 2 (II) gauge to vanishing quark mass; compare Fig. |7] 
with Fig. |^. In Fig. |?] we see large errors in the infrared 
region and points along the temporal axis lying below 
the bulk of the data. These are indicators of finite vol- 
ume effects, an unexpected result given that earlier gluon 
propagator studies ]9, 10 appear to conclude that Lapla- 
cian gauge is less sensitive to volume than Landau gauge. 

c? 2 (Ill) performs very poorly: see Figs. @ and 0. The 



FIG. 4: Mass function for quark mass ma = 0.05 (m ~ 76 
MeV), KS action (top) and Asqtad action (bottom) in Landau 
gauge. 



gauge fixing procedure failed for four of the configura- 
tions and eight of the remaining configurations produced 
Z and M functions with pathological negative values. We 
have seen that this type of gauge fixing fails to produce 
a gluon propagator that has the correct asymptotic be- 
havior jllj]. It seems likely that we are dealing with ma- 
trices with vanishing determinants, which are destroy- 
ing the projection onto SU(3). We expect the degree to 
which this problem occurs to be dependent on the sim- 
ulation parameters and the numerical precision used (in 
this work the gauge transformations were calculated in 
single precision). 



V. GAUGE DEPENDENCE 

Now we investigate the quark mass and Z functions 
in Landau, d 2 (l) and c? 2 (II) gauges. Fig. |Io| shows the 
Z function for the Asqtad action in Landau and <9 2 (I) 
gauges. They are in excellent agreement in the ultraviolet 
- as they ought - but differ significantly in the infrared. 
The Z function in the Laplacian gauges is more strongly 
infrared suppressed than in the Landau gauge. There 
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FIG. 5: Mass function for quark mass ma = 0.0125 (m ~ 19 
MeV), KS action (top) and Asqtad action (bottom) in Landau 
gauge. In both cases, lowering the quark mass increases the 
amount of noise, but the Asqtad action seems to be affected 
less than the KS action. Note that the y-axis for the bottom 
figure starts below zero. 



may be a small difference in Z(q 2 ) between <9 2 (I) and 
<9 2 (II) gauges. 

In all cases the quark Z function demonstrates little 
mass dependence. Deviation of Z from its asymptotic 
value of 1 is a sign of dynamical symmetry breaking, so 
we expect the infrared suppression to go away in the limit 
of an infinitely heavy quark. 



In Fig. 11 we show the Z 
function in Landau gauge for the lightest and the heav- 
iest quark masses in this study. The two are the same, 
to within errors, although if we look at the lowest mo- 
mentum data we see that the point for the low mass lies 
below the high mass one. Fig. fl2| shows Z in 9 2 (II) gauge 
for three quark masses. Again, the data are consistent, to 
within errors, but there is a systematic ordering of light- 
est to heaviest. We conclude from this that behavior is 
consistent with expectations, dependence on the quark 
mass - if any - is very weak. One possible explanation 
is that all the masses studied are light - less than or ap- 
proximately equal to the stange quark mass - and that 
heavier masses will affect the Z function more clearly. 



FIG. 6: Mass function for quark mass ma = 0.05 (m ~ 76 
MeV), Asqtad action in d 2 (I) and d 2 (II) gauges. Comparing 
with Fig. W we see that Landau, <9 2 (I) and <9 2 (II) gauges yield 
similar results. 



The mass functions in Landau and 9 2 (I) gauges - 
shown in Fig. [l3] - agree to within errors. The data for 
d 2 (l) gauge seems to sit a little higher than the Landau 
gauge through most of the momentum range, so with 
greater statistics we may resolve a small difference. The 
mass functions are nearly identical in d 2 (l) and c? 2 (II) 
gauges: see Fig. n4[ 



VI. MODELLING THE MASS FUNCTION 

The Asqtad quark mass function at each value of the 
mass has been cylinder cut and extrapolated - by a 
quadratic fit - to zero mass. The quadratic fit was cho- 
sen on purely practical grounds and a linear fit worked 
almost as well. A fit to each of the mass functions was 
then done, using the ansatz 



M(q) 



cA 1 



q 



2a 



A 2a 



m , 



(17) 



which is a generalization of the one used in Ref. 0. As 
we have seen, the quark mass function in the Laplacian 
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FIG. 7: Mass function for quark mass ma = 0.0125 (m ~ 19 
MeV), Asqtad action in c? 2 (II) gauge. We see a lot of infrared 
noise at this low quark mass in this gauge. Compare with the 
Landau gauge result in Fig. U. 



gauge is almost indistinguishable from that in the Lan- 
dau gauge, except that it is somewhat noisier (this may 
change at smaller lattice spacing). For this reason, we 
only show fits in the Landau gauge. Table I shows a 
sample of the fits. The table is divided into two sections, 
one in which the parameter a was held fixed, one in which 
it was allowed to vary. We see that for the heaviest mass, 
a = 1.0 provides an excellent fit, but in the chiral limit, 
a > 1.0 is somewhat favored. The role of a may be seen 
if Fig. |l5|. At a = 1.5, the infrared and ultraviolet are 
significantly flatter, while the mass generation at around 
one GeV is made steeper. 



TABLE I: Best-fit parameters for the ansatz, Eq. (|l7|), in Lan- 
dau gauge, in physical units. Where no errors are indicated, 
the parameter was fixed. 



m 


c 


A 


mo 


a 


M(0) 


2 

X 


(MeV) 




(MeV) 


(MeV) 




(MeV) 


per d.o.f. 


114 


0.35(1) 


910(20) 


142(7) 


1.0 


462(9) 


0.38 


95 


0.36(5) 


880(70) 


117(7) 


1.0 


440(20) 


0.42 


76 


0.39(5) 


830(70) 


92(7) 


1.0 


420(20) 


0.42 


57 


0.45(4) 


770(50) 


70(7) 


1.0 


410(20) 


0.51 


38 


0.49(8) 


720(60) 


44(6) 


1.0 


400(30) 


0.56 


19 


0.54(9) 


670(60) 


18(6) 


1.0 


380(30) 


0.69 
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FIG. 8: The quark Z function with the Asqtad action in 
c? 2 (Ill) gauge with ma = 0.05. Data represents 76 configu- 
rations (top figure), some of which actually provide negative 
contributions. The bottom figure shows the same data with 
the configurations producing negative contributions removed 
from the sample (8 configurations). 



In this model, a is acting as a function of the bare mass, 
controlling the dynamical symmetry breaking. Unfortu- 
nately, the paucity of data points in the infrared leaves a 
poorly determined in the chiral limit. Furthermore, the 
degree to which our infrared data may be affected by the 
finite volume and by the chiral extrapolation is not really 
known. Finally, this ansatz is still crude in that it does 
not provide the correct asymptotic behavior. 

As was explained is Section II, the quark mass function 
approaches the renormalised quark mass in the ultravi- 
olet, which itself becomes the bare mass in the q — ► oo 
limit. It is a general feature of this study that the ultra- 
violet tail of M{q) sits somewhat higher than the bare 
mass. The situation is summarised in Table |lj]. This de- 
viation from the correct asymptotic behavior is probably 
the consequence of an insufficiently small lattice spacing. 
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FIG. 9: The quark mass function with the Asqtad action in 
c? 2 (Ill) gauge with ma — 0.05. The upper data represents 76 
configurations, some of which actually provide negative con- 
tributions. Signal is almost completely lost. The lower data 
represents 68 configurations. Removing the negative contri- 
butions has only barely restored the signal. 



FIG. 10: Comparison of the quark Z functions for the Asqtad 
action for quark mass ma = 0.05. Points marked with open 
circles are in Landau gauge, solid triangles are in <9 2 (I) gauge 
and open triangles are in c? 2 (II) gauge. Data has been cylinder 
cut. 
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VII. CONCLUSIONS 

We have seen that the Asqtad action provides a quark 
propagtor with improved rotational symmetry compared 
to the standard Kogut-Susskind action and that the dif- 
ference between them increases as we go to lighter quark 
masses. Furthermore, we have seen that the Asqtad ac- 
tion has smaller mass renormalization and better asymp- 
totic behavior. 

Our results for the quark propagator show that the 
quark mass function is the same, to within statistics in 
d 2 (l) and 9 2 (II) gauges, while the Z function is slightly 
different. There is little difference between the quark 
mass function in Landau gauge and in Laplacian gauge, 
but the Z function dips more strongly in the infrared in 
Laplacian gauge than in Landau gauge. The infrared re- 
gion of the Laplacian gauge mass function seems to be 
particularly badly affected by decreasing quark mass. We 
have seen that the c? 2 (Ill) gauge gives very poor results 
in SU(3), in calculations of the quark propagator, consis- 



FIG. 11: Comparison of the Landau gauge quark Z functions 
for the heaviest (ma = 0.075) and the lightest (ma = 0.0125) 
quark masses, with the Asqtad action. Data has been cylinder 
cut. We see that over this range of values the mass depen- 
dence of the Z function is very weak. 



tent with results for the gluon propagator pJ| . Overall 
the Landau gauge results of this work for the Asqtad 
appear to be consistent within errors with the results 
of ealier Landau gauge studies of the quark propagator 
[4,5,6]. Our results suggest that the M function is in- 
sensitive to whether we use Landau or Laplacian gauge, 
whereas the Z function has an enhanced infrared dip in 
Laplacian gauge. 

As we have simulated on only one lattice, it remains 
to do a thorough examination of discretization and fi- 
nite volume effects. The chiral limit was obtained by 
extrapolation, which may provide another source of sys- 
tematic error. It will also be interesting to investigate the 
0(a 4 ) errors by using an improved gluon action. Further 
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FIG. 12: Comparison of the quark Z functions for the three 
quark masses ma — 0.0125, 0.025 and 0.05, with the Asqtad 
action in <9 2 (II) gauge. Data has been cylinder cut. As in 
Landau gauge, they agree to within errors, although there 
is a systematic ordering of the infrared points from heaviest 
quark (top) to lightest (bottom). 
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FIG. 14: Comparison of the quark mass functions for the 
Asqtad action, ma = 0.05. Points marked with solid triangles 
are in d 2 (V) gauge and open triangles are in c? 2 (II) gauge. Data 
has been cylinder cut. 



TABLE II: Estimates of the Landau gauge quark mass func- 
tion at zero four-momentum. Here m is the bare input mass 
and mo is the ultraviolet mass from the fit of the mass func- 
tion. 





a = 


1.0 


a ~ 1.3 




m (MeV) 


M(0) (MeV) 


mo/m 


M(0) (MeV) 


mo/m 


114 


462(9) 


1.25(6) 


428(7) 


1.35(6) 


95 


440(20) 
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76 
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FIG. 13: Comparison of the quark mass functions for the 
Asqtad action, ma = 0.05. Points marked with open circles 

are in Landau gauge and solid triangles are in <9 2 (I) gauge. APPENDIX A- STAGGERED QUARK 

Data has been cylinder cut. PROPAGATORS 



work may allow the development of a more sophisticated 
ansatz which has the correct asymptotic behavior. Fi- 
nally, studies on finer lattices could be used to calculate 
the chiral condensate, light quark masses and potentially 
Aqcd- 
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In this appendix we give details of the quark propa- 
gator calculation using the Kogut-Susskind and Asqtad 
actions. The free KS action is 

+ m^xWxW ( A1 ) 

X 

where the staggered phases are: tj^x) = (— 1)^ A ' x and 

c M = (lXv<H (A2) 
L otherwise. 
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where the 7„ satisfy 



0.30 
0.25 

o.ao 

0.15 

^ o.io - 

0.05 - 
0.00 

-0.05 





q (GeV) 



FIG. 15: Mass function extrapolated to the chiral limit. Er- 
rors are Jack-knife. Fit parameters are top: c = 0.080(20), 
A = 520(50) MeV, m = 0.0, a = 1.0, x 2 I dof = 1.3 and 
bottom: c = 0.030(4), A = 870(60) MeV, m = 0.0, a = 
1.52(23), x 2 I dof = 0.49. 



To Fourier transform, write 
27m„, 



as fc M = p p + ira^ , where 

27TTO,, 



n„ = 0, . . . , L M - 1 (A3) 



L„ 



m„ = 0,...,^-l 



= 0, 1, 



and define f, = i ^Zfc- Then 



E 



X (a0= / e*- X (fc)= / E e ' 



(p+7ra)-a: 



Xa(p)- 



Defining 



<5q/3 — II^^^I m od 2 

(7^ = (-1)"" 



(A4) 
(A5) 



(A6) 
(A7) 



(A8) 
(A9) 



+ T 1 * 



(A10) 
(AH) 



forming a "staggered" Dirac algebra. Putting all this 
together, we can derive a momentum space expression 
for the KS action. 



J P a/3 ll 



(A12) 

From this we can see that, in momentum space, the 
tadpole improved, tree-level form of the quark propagator 
is 

S~l& m) = u i (7 M ) Q/3 sin(p M ) + mS a/3 (A13) 



wher e p u is the discrete lattice momentum given by 
Eq. ( |A4| ) . Assuming that the full propagator retains this 
form (in analogy to the continuum case) we write 

S ^(P) = * E fiJap s ™(pMp) + B{p)~5 a p (A14) 
= Z-\ P ) [i (7 M ) Q/3 sin(p M ) + M (p)S a 



For the KS action, it is convenient to define 

= sin(p M ) 



■ 
(A15) 

(A16) 



as a shorthand. 

Numerically, we calculate the quark propagator in co- 
ordinate space, 

G{x,y) = (x{x)x{y)) 

= / exp{i(p + ira)x — i(r + 7r/3)y} 

a/3 JP' r 

X (Xa(p)Xp(r)) 
— / exp{i(p + ira)x — i(r + 7r/3)y} 

a/3 J* > r 

X 5 pr Sa/3(p) 

= E / exp-^K^ -y)} 



i/3 •'P 



x exp{i7r(ax - /3y)}S afj (p). (A17) 

To obtain the quark propagator in momentum space, we 
take the Fourier transform of G(x, 0) and, decomposing 
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the momenta as = + 71 5^ 1 < r M < 7r, get 
G(fc) = G(r + tt<5) = G«(r) = e- lfcE G(izi, 0) 

a/3 J P x 

x exp{i(p + na)x}S a p(p) 

= Z2 SprSaSSapip). (A18) 
a/3 "'P 

Thus, in terms of the KS momenta, q, 



from which we obtain 



^TrG a (g) = 16JV, 



q 2 + M 2 (q) 



Z(q)M(q) 



(A19) 



q 2 + M 2 (q) 
lW c B(q), (A20) 



and 



i££(-l)N„Tr[G 5 ( 9 )] = 16N ^ 2 q ^M2 iq) 



s P 

Note: we could determine 
A{p) = Z-\p) 



l6N c q 2 A(q). (A21) 



B(p) = 



16iV c E,sin 2 (^) 

x EE(%) a ^ sin fe) Tr K(p)] ( A22 ) 

a/3 p 

M(p) 1 



167V, 



(A23) 



but we would rather avoid inverting S a p(p). 
Putting it all together we get 



A(q) = Z-\q) = 
B(q) = 
M(q) = 



A(q) 



A 2 {q)q 2 +B 2 (q») 
M(q) _ B(p) 
Z(q) A 2 (q)q 2 + B 2 (p) 
B(q) 



A( q y 



(A24) 
(A25) 
(A26) 



The tadpole improved, tree-level behavior of the Z and 
mass functions are simply 



u 



(A27) 



and 



M u = 



(A28) 
respectively. 

The Asqtad quark action S is a fat-link Staggered 
action using three-link, five-link and seven-link staples 
along with the three- link Naik term |l3) and five-link Lep- 
age term H], with tadpole improved coefficients tuned 
to remove all tree- level 0(a 2 ) errors. This action was mo- 
tivated by the desire to minimize quark flavour changing 
interactions, but has also been reported to have good 
rotational symmetry. 

At tree-level (i.e. no interations, links set to the iden- 
tity), the staples in this action make no contribution, so 
the action reduces to the Naik action. 



S = \^x{x)r)n{x) g (x{x + n) - x(x - m)) 



1 

24 



(A29) 



The quark propagator with this action has the tree-level 
form 



S a p(p; m) = u Q i £ (7 M ) Q/3 sin(p M ) - sin(p M ) 



- — sin(3p M ) + m8 a p (A30) 



so we choose 



9 1 

«V»Gv) = g sin(p M ) - — sin(3p M ) 

= sin(p M ) 1 + isin 2 (p AI ) 



(A31) 
(A32) 



Having identified the correct momentum for this action, 
we can calculate the invariant functions as before. No 
further tree-level correction is required. 
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